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In the next decade or two, extremely accurate tests of general relativity under extreme con- 
. ditions are expected from gravitational wave observations of binary black hole mergers with a 

wide range of mass ratios. In addition, major improvements are planned in both strong and 
weak equivalence principle tests; clock measurements based on the ACES program on the ISS; 
more accurate light-bending measurements; and other new types of tests. However, whether 
these tests are all consistent with general relativity or not, it still appears desirable to proceed 
with a much improved measurement of the Shapiro time delay. A suggested approacff is 
based on using a high-quality optical clock in a drag- free spacecraft near the sun-earth LI 
point and a smaller drag- free transponder spacecraft in a two-year period solar orbit. Laser 
phase travel-time measurements would be made between the two spacecraft over a period of 10 
or 20 days around the time when the line of sight passes through the Sun. The requirements 
on the optical clock stability and on the drag- free systems will be discussed. The accuracy 
achievable for the time-delay appears to be better than 1 part in 100 million. 
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1 Introduction 

- '— j . The first suggestion to measure the gravitational time delay for electromagnetic waves passing 

r> \ near the sun was made by Irwin I. Shapiro in 1964^ The extra gravitational time delay for 

' two-way measurements of light propagating from Earth to a spacecraft passing behind the Sun 

can be more than 200 microseconds. In the Parametrized Post-Newtonian (PPN) formulation of 
gravitational theory, the main contribution to the time delay is proportional to (1 -|- 7), where 
7 is a measure of the curvature of space. In General Relativity (GR), 7 = 1. 

In view of the well-known lack of a theory that connects GR with quantum theory, im- 
provement of high-accuracy tests of the predictions of GR should be the object of research in 
the coming decade. Many alternatives to GR involve additional scalar fields. Studies of the 
evolution of scalar fields in the matter-dominated era of the universe indicate that the universe's 
expansion tends to drive the scalar fields toward a state in whic h the scalar-tensor theory is only 
slightly different from GR. Some scalar-tensor extensions of GRpH! predict deviations from the 
GR value of 7 in the range from 10~^ to 10~®. Improved information about 7 would provide 
important insight into the evolution of the universe and directly limit the range of applicabilituy 
of alternative gravitational theories. 

Recently, a measurement of 7 with accuracy ±2.3 x 10~^ was made during the Cassini 
missionP Further improvements in the accuracy for 7 to roughly 10~^ are expected from two 
missions of the European Space Agency (ESA): the GAIA astrometric mission, which will mea- 
sure the gravitational deflection of light rays by the sun, and the Bepi Colombo mission to 



Mercu ry, wh ich will make improved measurements of the solar time delay. We describe here a 
missioiP^^ that can reach an accuracy of about 1 x 10~* for determining 7. 

2 Mission orbits and predicted time delay 

For the proposed mission, one spacecraft (SI) containing a highly stable optical clock would be 
placed in an orbit near the LI point, about 1.5 million km from the Earth in the direction of 
the Sun. The second spacecraft (S2) would have a 2 year period orbit in the ecliptic plane, with 
an eccentricity of 0.37. S2 would pass through superior solar conjunction about 1, 3, and 5 years 
after launch and would be near aphelion at those times. Both spacecraft would have drag-free 
systems to nearly eliminate the effects of spurious non-gravitational forces. A measurement of 7 
to a level of 1 x 10~^ would be carried out by observing the time delay of laser signals exchanged 
between the two spacecraft when the line of sight passes near the Sun's limb. Atmospheric 
effects would be absent and continuous observation would be possible. With S2 near aphelion, 
the range rate would be low, and the orbit determination problem would be much reduced. 

The crucial measurements of time delay occur within a few days of superior conjunction 
and are primarily characterized by a logarithmic dependence on the distance of closest approach 
of the light to the mass source. The predicted gravitational time delay due to a non-rotating 
mass source, expressed in terms of the radii r^, rs of the endpoints of the photon path, and the 
elongation angle $ between the radius vectors from the source to the endpoints, is^^ 

^ ^l frA + rB+rAB\ + 7?rAB , Aab(8 - 4/3 + 87 + 3e) 

cAtdelay = U + 7 log ■ — — + — (1 

where // = GMq/c^, and tab is the geometric distance between the endpoints in isotropic 
coordinates and /3 and e are PPN parameters measuring the nonlinearity of the time-time and 
space-space components of the metric tensor. In GR, (8 — 4/3 -|- 87 -|- 3e)/4 = 15/4. The time 
delay in Eq. ([T]) is expressed in terms of observable quantities, and does not involve the unknown 
impact parameter or distance of closest approach. The non-linear terms are a few nanoseconds 
so they are significant, but do not have to be estimated with great accuracy. The contributions 
to time delay due to the solar quadrupole moment are small and can be estimated with sufficient 
accuracy that they will not contribute significantly to the error budget. 

The measurements will be made by transmitting a laser beam with roughly 40 GHz sidebands 
on it from from SI to S2, and comparing with a similar beam generated on S2 and sent back 
to SI. From the phase differences of the sideband beat notes, the round-trip delay time can be 
obtained. With 20 cm diameter telescopes, and given the one-way travel time of about 1600 
s, the received signal would be roughly 1000 counts/s for 1 W of transmitted power. This is a 
weak signal, but it is strong enough so that the chances of a cycle slip should be very small. If 
we consider the round-trip delay times Atdeiay to be the observable, then the change in delay 
from 0.75 days to 4 days on either side of conjunction is about 64 microseconds. 

3 Signal-to-Noise Analysis 

We can estimate the lowest possible uncertainty that could be attained in this experiment on 
the basis of the optimal Wiener filter, which t akes a dvantage of the known time signature of 
the signal and includes the expected clock noiseP^^ For this case, uncertainties in the various 
orbit parameters for the two spacecraft are ignored, and the travel time between the spacecraft 
is assumed to be constant except for changes in the gravitational time delay. The time signature 
of 7* = (1 -|- 7)/2 is taken to be represented by the logarithmic function 



git) = -B{log\Rt\-M) 
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where M is the mean value of log \ Rt\ over the time periods —t2 to —ti and ti to t2 (a short time 
interval during occultation is excluded), and for the proposed experiment B = 0.97 x 8/i/c = 
3.82 X 10~^ s. The rate at which the line of sight to the distant spacecraft passes across the sun 
is i? = 1.9 solar radii per day. 

Let g{f) be the Fourier transf orm o f g{t) over the time of the measurements. Then the 
signal-to-noise ratio may be foundU™ in terms of an integral over all frequencies of 
An important consequence of the logarithmic form of the time delay, Eq. (2), is that if the 
noise has a constant spectral density, only about 2.5% of the signal-to-noise ratio comes from 
frequencies below 1 microHz, where the acceleration and clock noise are expected to increase. 
Just integrating down to 1 microHz, we find an uncertainty less than 1 x 10~^ for 7. Almost all 
of the power in \g{f)\'^ is at frequencies between 1 and 8 microHz, so it is clear that the noise 
at these low frequencies will provide the main limitation on the results. 

Actually, all of the in-plane parameters for the orbits of the two spacecraft have to be 
solved for, as well as 7. Our model for this includes uncorrelated 0.02 picosecond uncertainties 
for measurements of the round-trip travel time over 3-hr periods. This is in addition to our 
assumed white clock frequency noise of 5 x 10~^^/\/Hz down to at least 1 microHz. Spurious 
acceleration noise is not included in this model, but its effect has been estimated to be small 
with our assumptions about its spectrum. The resulting uncertainty in 7 is less than 1 x 10^^. 

So far, we have assumed that time-delay measurements are only made over a total period of 
8 days around solar conjunction. This was done in order to make sure that spurious acceleration 
noise at frequencies below 1 microHz would have little effect. However, simulations for longer 
observing times are desirable, with full allowance for spurious acceleration noise at the lowest 
frequencies, as well as for the orbit determination part of the problem. The longer observation 
period may help to improve the determination of orbit parameters as well as 7. 

4 Spacecraft SI clock 

The major requirement for the mission is to fly an optical clock on SI that has very high 
stability over a period of at least 8 days around superior conjunction. The nominal design goal 
for the mission is to achieve a fractional frequency noise power spectral density amplitude of 

5 X 10^^^/\/Hz from 1 Hz down to at least 1 microHz. (This is nearly equivalent to an Allan 
deviation of 5 x 10^^^ / ^/t for times from 1 s up to 10® s.) 

As an example of the desired performance, a spectral amplitude of about 2 x 10~^^/\/Hz has 
been achieved in the laboratory down to 1 mHz for the 267 nm transition in sympathetically- 
cooled A1+ ions inamagnetic trapff^ Other l eading candidates for optical clocks in space are 
cooled Sr^s atomPll^ and Yb^^i atomPISl 

in optical latices. However, substantial develop- 
ment is needed to show that such optical clocks can be designed for use in space and can be 
space qualified. 

5 Drag-free system 

The required performance builds on that planned for the LISA mission. For frequencies down to 
10~^ Hz for LISA, the requirement on the acceleration power spectral density amplitude is less 
than 3x 10~^^m/s^/\/Hz. However, the performance is expected to degrade at lower frequencies. 
The main challenge for achieving good performance at low frequencies is minimizing thermal 
changes, and particularly thermal gradient changes, near the freely floating test mass in the 
drag-free system. On LISA this is done almost completely by passive thermal isolation. For a 
time delay mission, a fairly slow active temperature control system would be used at frequencies 
below 10~^ Hz. Changes in solar heat input over the 8 days around conjunction would be quite 



small for S2, because conjunction occurs near aphelion. The required drag-free performance is 
roughly 1 x 10~^^m/s^/\/Hz down to 1 microHz. 

In fact, much of the desired freedom from spurious accelerations needed for LISA has been 
demonstrated in the laboratory with torsion pendulum measurements?^ But, more important, 
the overall performance of the drag-free system will be dem onstrated in the LISA Pathfinder 
Mission, which is scheduled for launch by ESA in 2014™ 

6 Other scientific benefits from the mission 

Additional effects such as those arising from non-linear terms in the 00-component of the metric 
tensor, parameterized by /?, as well as other time delay effects originating in the sun's rotation, 
can also be measured. The clock at the LI point will experience frequency shifts from the earth's 
potential, solar tidal effects, and second-order Doppler shifts. Relative to a reference on earth's 
surface, the fractional frequency shift is about -1-6.9 x 10^^*^, and is almost all gravitational. 
Comparing the clock at LI with a similar clock on earth's geoid will give accuracies of a few 
parts per million in a few hours, which is orders of magnitude more accurate than the Vessot- 
Levine 1976 Gravity Probe A result. This result is comparable to that expected from the 
upcoming ACES mission?^ 

7 Postscript 

After the Moriond Meeting we learned about proposalP^^^Q] for a mission called ASTROD I, 
with improved measurement of the gravitational time delay as one of its main objectives. In the 
proposed mission, the time delay measurements would be made between a drag-free spacecraft 
in a solar orbit with a semi-major axis of about 0.6 AU and laser ranging stations on the Earth. 
The projected accuracy for determining the PPN parameter gamma is 3 x 10~^. 

In these papers the drag-free requirement is given explicitly only over the frequency range 
from 0.0001 to 0.1 Hz, and is 3 x 10^^^ m/s2/\/Hz at 0.0001 Hz. And the only clock frequency 
stability requirements given are 1 x 10~^^ for the clock in the satellite by comparison with 
ground clocks and 6 x 10~^^ stability in the round-trip travel time of roughly 1700 s. However, 
to reach the accuracy goal given for gamma appears to require low levels of spurious acceleration 
noise and clock noise down to about 1 microHz or lower. Thus it seems possible that quite low 
spurious acceleration and clock noise levels at low frequencies actually were implemented in the 
simulations on which the ASTROD I accuracy goals are based. 
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